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Breaking down the entire W boson spin observables from its decay 
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We discuss the eight independent spin observables for the W boson in terms of its vector and 
tensor polarizations and identify the angular distributions and asymmetries able to separate them in 
collider experiments. The results are applied to the study of polarized top quark decays and diboson 
resonances. These novel observables are of great value for disentangling new physics mechanisms in 
W boson production. 


The discovery by the ATLAS [l| and CMS 0 experi¬ 
ments at the Large Hadron Collider (LHC) of a new par¬ 
ticle with a mass of 125 GeV and the expected properties 
of the Higgs boson has marked the programme of high 
energy physics for the next coming years. On one side, 
it is mandatory to be precise enough in the measurement 
of its properties in order to ascertain its nature as rem¬ 
nant of the electroweak symmetry breaking mechanism. 
On the other hand, the search of new physics beyond the 
standard model (SM), related or not to this symmetry 
breaking problem, makes the W boson —as mediator of 
the weak charged interactions— a privileged particle to 
be a messenger of the scenario chosen by Nature. In ad¬ 
dition, new unknown production mechanisms will leave 
their imprint in its detailed properties. 

Interestingly enough, being a spin-1 particle, the W 
boson provides 8 additional observables, beyond the 
number of events, at our disposal to probe new physics 
and even discriminate among new physics models. As a 
consequence, it is of the highest importance to give an 
unambiguous procedure for extracting these spin proper¬ 
ties from the experiments. This is the main motivation 
of this work. Some partial treatments have appeared in 
the literature 0, and also for specific reactions, like e+e“ 
annihilation 0 , but a direct methodology for connecting 
the decay distributions to the three polarizations and five 
alignments —or tensor polarizations— of the W boson is 
lacking. Our results will be first applied to the decay 
of the top quark, which is the heaviest known particle 
and thus widely regarded as a promising probe of new 
physics im. Subsequently, we will consider the pro¬ 
duction of a WZjWW diboson resonance, motivated by 
a recent 3.4cr excess reported by the ATLAS Collabora¬ 
tion 0. 

The role of the various mechanisms for the produc¬ 
tion of the W boson is to describe its spin state in the 
form of a 3 X 3 density matrix, Hermitian with unit trace 
and positive semidefinite. By fixing a coordinate system 
(x, y, z) in its rest reference frame, it can be written in 
terms of the three spin operator components Sm and the 


five tensor operator components Tm as 
1 2 

P=l'^ + \ X {SmYSm+ X (TaiYTm, (I) 

M=-l M=-2 

with S±i = T^-^{Si±iS2), So = S3 the spin operators in 
the spherical basis and Tm five rank 2 irreducible tensors, 
which are built from Sm as 

t ±2 = SI 3 , r±i = ^ [s±iSo + SoS±i], 

To = ^ [S+iS-i + S-iS+i + 2SY\ . (2) 

Their expectation values {Sm) and (Tm) are the multi¬ 
pole parameters corresponding to the three polarization 
and hve alignment components. More conveniently, the 
second term in Eq. © can be rewritten using the spin op¬ 
erators in the Cartesian basis, and the third one defining 
the Hermitian operators 

Ai = i(Ti - T_i), A2 = i(Ti + T_i), 

Hi = i(T2+T_2), H2 = i(T2-T_2). (3) 

Explicitly, the W boson density matrix parametrized in 
terms of expected values of observables reads 

P±i±i = I ± ^(Sa) + -^{To) , 

P±w = ^ [(-5i) T i{S2)\ T [(Ai) T *(^ 2 )] , 

Poo = \ - -^(To ), Pi -1 = (Hi) - i(H2), (4) 

and Pm'm = Ptnm'- This density matrix entirely deter¬ 
mines the angular distribution of the W decay products 
in —i £^iy, with i = e,p. Using the helicity formal¬ 
ism of Jacob and Wick 0, the amplitude for the decay 
of a IT boson with third spin component m into charged 
leptons £, v with respective helicities Ai, A 2 is written as 

MmAiA. =aA,A.H),:;,((/.*,r,0), (5) 

with {0*,(j>*) the polar and azimuthal angles of the 
charged lepton momentum in the W boson rest frame, 
A = Ai - A 2 and H^/^(a,/3,7) = e"“'"'e-*T'™(i^,^(/3) 


2 


the so-called Wigner D functions ; are con¬ 

stants, and all the dependence of the amplitude on angu¬ 
lar variables of the final state products is given by the D 
functions. For the SM left-handed interaction of the W 
boson with the charged leptons, and taking them mass¬ 
less, one has (Ai,A 2 ) = (± 1 / 2 ,=f 1 / 2 ) for —>■ 

Hence, the fully differential W decay width can be writ¬ 
ten as 


dr 

dcos9*d(j)* 


m,m' 


( 6 ) 

with C a common prefactor including the phase space 
and the single non-zero Here, A = 1 (—1) for 

bF+ {W~) decays and the functions d^^,^{l3) are real, 
and so is the right hand side of Eq. Explicitly, the 
normalised distribution reads 


f = I? {5(1 + <=“"«•) + 


{T«) 


+ {Si) cos (j)* sin 9* + {S 2 ) sin 0* sin 9* 

— (Hi) cos (j)* sin 29* — {A 2 ) sin (p* sin 29* 

. / R. \ rtr\ci O / R„\ oin OAt* 


-f (Hi) cos 2 ^ 1 * 


(7) 


An analogous expansion has been used by the CDF [TJ 
and CMS [l2| Collaborations to study the dilepton an¬ 
gular distribution in Z production, but an identification 
of the coefficients of the different terms with spin observ¬ 
ables has not been provided — indeed for the Z boson 
the relation is more contrived since both dilepton he- 
licity combinations (Ai,A 2 ) = (± 1 / 2 ,=f 1 / 2 ) contribute 
to the decay. This expansion is also implicit in Q but 
the correspondence with W boson spin observables is not 
made (instead, predictions for some asymmetries are ob¬ 
tained for the specific case of e+e“ —>■ W'^W~ with gen¬ 
eral couplings). Polar asymmetries have been proposed 
in 0, to measure (Hi) and (H 2 ). Azimuthal asymmetries 
have also been proposed to test the spin of the decay¬ 
ing particles [T^ [ij, in particular the spin -1 nature of 
W and Z bosons can be tested by measuring the sin 2(j)* 
dependence M- (Within our framework, this amounts 
to measuring non-zero (Hi) or (H 2 ), which clearly shows 
that the spin is larger or equal than one, given the rank of 
the irreducible tensor T 2 .) But, as we have emphasized, a 
general method for connecting all W boson polarization 
observables with angular distributions and asymemtries 
has not yet been provided. 

The distribution © can in principle be used to de¬ 
termine all the W boson spin observables with an eight- 
parameter fit. However, it is much simpler to project on 
certain directions to obtain one-dimensional distributions 
that depend on fewer parameters. This also allows for a 
cleaner extraction from experimental data without any 
loss of information. We do this in detail for the decay 
of a W~^ boson, and point the differences for W~ boson 


decays. The asymmetries obtained here could as well be 
applied to the determination of all the parameters in the 
Z boson distribution, extending current analyses [l^ to 
the full set of eight coefficients. 

Integrating over the azimuthal angle, the off-diagonal 
contributions vanish and the polar angle distribution 
reads 


1 dr 
F dcos 9* 


|h+ (1 + cos 0*)2 -t ^Ho sin^ 9* 

8 4 

+ ^F_{1-cos9*)\ (8) 


where F+ = pn, Hq = poo and H_ = p_i_i are the well- 
known W polarisation fractions [3, [3 1 which are not in¬ 
dependent but satisfy F+ -f Hq -I- H_ = 1. In the common 
case that the z axis is taken in the direction of the W mo¬ 
mentum, F±fi are the helicity fractions. These have been 
measured by several experiments for the W bosons pro¬ 
duced in top quark decays [TtI - I^ . But instead of mea¬ 
suring three quantities that are not independent, it would 
be simpler to just use two independent ones: the longi¬ 
tudinal polarization (S 3 ) and the alignment (Tq), which 
obviously contain the very same information. Writing 
the explicit values of pn, poo and p_i_i from Eqs. ®, 
the resulting distribution is 


1 dr 

F dcos 9* 


1 + 3 (^ 3 ) cos r + ^(ro) [3 cos2 0 *- 1 ] . 


(9) 

With a two-parameter fit to the cos 9* distribution, (S 3 ) 
and (To) can be measured as well. For illustration, 
the functions multiplying these observables are shown in 
Fig. [ 1 ] (left). Alternatively, the two observables may be 
determined from asymmetries in the cos 9* distribution: 
the forward-backward (FB) asymmetry Hfb [H, [ 2 ^ and 
an “edge-central” asymmetry Hec, 


Hfb = i[r(cos0* >O)-r(cos0* <0)] = ^(Sa), 
Hec = i[r(|cos0*|>i)-r(|cos0*|<i)] 

= ( 10 ) 


Whether this simpler approach also gives more precise 
results from the experimental point of view, once the 
experimental uncertainties are taken into account, has to 
be investigated with dedicated measurements. 

In order to measure off-diagonal density matrix ele¬ 
ments one must replace the integration over (p* by an 
integration using a measure f{p*) that avoids the can¬ 
cellation in the total integral, while selecting the de¬ 
sired observables. We first use fi{(p*) = sign cos ())*, 
f2{<P*) — sign sin ^*. Experimentally, integrating with 
fi amounts to introducing a weight of -fl for events 
with (p* G ] — 7 r/ 2 , 7 r/ 2 [, and of —I for events with 
p* G ] 7 r/ 2 , 37 r/ 2 [. Integrating with /2 amounts to intro¬ 
ducing a weight of -1-1 for events with p* G ]0, 7 r[, and of 
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FIG. 1: Left: functions multiplying {S 3 ) and (To) in the distri¬ 
bution and constant term. Right: functions multiplying 
{Si, 2 ) and {Ai, 2 ) in the distributions (Hill . 


— 1 for events with (p* G Jtt, 2tt[. The resulting quantities, 
which we denote as iJiT, (52r, have an angular distribution 

with k = 1,2. The functions multiplying {Sk) and {Ak) 
are shown in Fig. [T] (right). These four observables can 
be measured by fitting the above distributions or, in a 
simpler way, by either integrating over cos 9* or doing 
the forward-backward asymmetry in this angle. The total 
integrals on the polar angle 6* give 

^ = 1 [T{cosr > 0) - T{cosP* < 0)] = ^{Si ), 

^ ^ [r(sin())* > 0) - r{sm(j)* < 0)] = ^( 52 ). 

( 12 ) 

Notice that (JiT/T and (52r/r are merely the FB asym¬ 
metries with respect to the x and y axes, respectively Q, 
so the relation with {Si) and {S 2 ) is expected. The FB 
asymmetries Apg in the ^fcF distributions m yield the 
other two observables, 

^FB = ^ [r(cOS(/)* C0S6>* > 0) — r(cOS()>* COS0* < 0)] 

= --(^ 1 ), 

TT 

^FB = ^ [r(sin(/)* COS0* > 0) — r(sin(/)* cos6** < 0)] 

= --(^ 2 ). (13) 

TT 


be performed to get azimuthal asymmetries 
A\ = [r(cos2(/)* > 0) - r(cos2(^* < 0)] = -{Bi ), 

= i [r(sin2<(.* > 0) - r(sin2r < 0)] = -{B 2 ). 

^ i TT 

(14) 


For W~ boson decays the differential distributions are 
obtained from the ones presented by replacing sin 6 >* —>■ 
— sin0*, COS0* —>■ — COS0*. The relations between asym¬ 
metries and observables in Eqs. m, m, o and m 
are the same for W decays, with extra minus signs for 
Afb and iJ^F/F, i.e. those asymmetries proportional to 
the expected values of the spin operators. 

Application to polarized top quark decays. The density 
matrix elements for W bosons resulting from the decay 
of a fully polarized top quark were given in Ref. Q for 
a general Wtb vertex including left and right vector and 
tensor couplings, using the helicity basis. Setting our z 
axis in the direction of the W boson momentum in the top 
quark rest frame, and the top spin direction in the x — z 
plane, s) = (sin 0 , 0 , cos 0 ), one has the (unnormalized) 
density matrix elements 


p±i±i = .So(l ± COS0) ± 2—i3i(l ± COS0). 

mt 


P±w — 


(Co + iBii) ± —(Cl + iDi) 


_y 2 Mi 


w 


\/ 2 M] 


w 


Poo — -4o + 2— Ai COS 9 , pi_i — 0, 
mt 


sin 0 , 
(15) 


and Prn'm = Pmm'^ with m* and Mw the top quark and 
W boson masses, respectively, g the modulus of the W 
boson 3-momentum in the top quark rest frame, and 
-4o,ii -Bop, Cop, T’op constants that in general depend on 
the Wtb couplings. Their explicit expressions were given 
in Ref. Q . (For top anti-quark decays one has to reverse 
the sign of the ,Bi, Vq, Ci and Ai terms.) For simplicity 
we consider here the quantities integrated over 9, despite 
an analysis dependent on this angle, or with 9 restricted 
to some interval, is possible with enough statistics. Inte¬ 
grating over cos 9 and normalizing to the W ^ £iy width, 
one has 


(- 53 ) = 4 


q 




isi) =; 


mt Ao + 2i3o ’ 
Co 


(^ 1 ) = - 


2 Mw Ao 

TT q 


-2Bo ’ 

Cl 


irp \ _ /2 Bq - Aq 

^ ~ i 3Ao + 2Bo 


4 Mw Aq 2Bo 


{S 2 ) = -t; 
{A2) = 


TT q 


Vi 


2 Mw Ao + 2Bo 

TT mt Vq 
4 Mw Aq -(- 2Bq 
(16) 


Let us now consider the measurement of (Bk). Inte¬ 
grating the doubly differential distribution with 51 (</>*) = 
sign cos 2(/)*, 32 (</>*) = sign sin 20*, all terms in Eq. (| 6 ]) 
vanish except those involving p±i=pi. Furthermore, 31 
selects the terms proportional to {Bi) and 32 the ones 
proportional to {B 2 ), so the integration over cos 0 * can 


and {Bi) = {B 2 ) = 0 since 3 ±ipi = 0 . 

High energy collider processes do not produce fully po¬ 
larized top quarks but with a partial polarization. For 
example, in t-channel single top quark production at 
the LHC, top quarks have a longitudinal polarization 
Pz ~ 0.9 in the direction of the spectator jet. In this 
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process, the top transverse and normal polarizations 
vanish [^, so one can take into account in Eqs. (ITSl) the 
degree of top quark polarization by simply introducing 
a Pz factor in the expressions of { 81 ^ 2 ) and (^ 1 , 2 )- We 
remark that P^ can be measured in a mo del-independent 
fashion [ 2 ^, so there is no need to assume that it takes 
the SM value in order to measure the W boson spin ob¬ 
servables. Four of these observables are nonzero in the 
SM, and their tree-level values for P^ = 1 are collected 
in Table m 



(5i> 

(^2) 

{S 3 ) {To) 

(^ 1 ) 

{A 2 ) 

SM 

0.510 

0 

-0.302 -0.445 

0.255 

0 

gn = 0.03 

0.500 

0 

-0.278 -0.472 

0.249 

0 

qr = 0.10 i 

0.507 

-0.084 

-0.284 -0.434 

0.253 

-0.042 


TABLE I: SM tree-level predictions for W spin observables 
in polarized top quark decays, and some predictions in the 
presence of an anomalous Wtb dipole interaction. 


The experimental values of ( 83 ) and (Tq) can be deter¬ 
mined from current measurements by using their relation 
with the W helicity fractions, 

(, 53 ) = F+-F., {To) = ^ [1 - 3Fo] . (17) 

Using the most precise measurements, obtained by the 
CMS Collaboration in ti production at 7 TeV [ 2 ^, we 
obtain { 83 ) = -0.302 ± 0.036, (Tq) = -0.427 ± 0.054, in 
good agreement with the SM prediction. (Next-to-next- 
to-leading order predictions for 0 are available [^ . 
but the differences with the tree level results are below 
the 0.3 ct level and for consistency with the rest of ob¬ 
servables we use the tree-level predictions in order to 
calculate { 83 ) and {To).) Notice that the uncertainty 
on { 83 ), obtained here by error propagation, would be 
around a factor of two smaller with a dedicated measure¬ 
ment. { 81 ) and (Ai), both nonzero in the SM, have not 
been measured so far. In this respect, we point out that 
the measurement of (^i) in single top quark production 
would be one of the very few, if not the first, measure¬ 
ment of the M > 0 tensor polarization of an elementarty 
particle.^ Remarkably, { 81 ), {Ai) are measurable with 
current LHC data at 8 TeV, using the same signal iden¬ 
tification and back grou nd rejection strategies of current 
measurements , but considering different angu¬ 

lar distributions. The expected uncertainties are around 
0.03 — 0.05, i.e. of the same order as those for { 83 ) and 


^ The OPAL and Delphi Collaborations measured the W 
density matrix elements at the Large Electron Positron collider, 
as a tool to search for anomalous triple boson couplings. The 
measurements of Repi_i and Impi_i could be intepreted as 
measurements of tensor polarizations but the uncertainties were 
so large that it is difficult to claim that a measurement was made. 


Resonance 

(Si) 

(S 2 ) (S 3 ) 

(To) 

(^ 1 ) (^ 2 ) 

(Bi) 

{B 2 ) 

WZ J = 0 

0 

0 

0 

-0.815 

0 

0 

0 

0 

WZ J = 1 

0.117 

0 

0 

-0.801 

0 

0 

0.003 

0 

WW J = 0 

0 

0 

0 

-0.815 

0 

0 

0 

0 

WW J = 1 

0.099 

0 

0 

-0.800 

0 

0 

0.003 

0 


TABLE II: Spin observables for bosons resulting from 
the decay of a heavy scalar or vector WZjWW resonance. 
Entries numerically smaller than 10~® are indicated as zero. 


(To). { 82 ) and {A 2 ) are nonzero only if there are anoma¬ 
lous imaginary contributions to the Wtb vertex; among 
them, only { 82 ) has been experimentally probed in t- 
channel single top quark production at 7 TeV by the 
ATLAS Collaboration [ 2 ^. Assuming the SM value for 
the top quark polarization (since it has not been yet mea¬ 
sured independently), the measurement of i52r/r yields 
{ 82 ) = —0.045 ± 0.106. Here the uncertainty is domi¬ 
nated by the statistical error, so a large improvement, by 
a factor of two or three, is foreseen with already collected 
8 TeV data. 

These polarization observables are also sensitive to new 
physics in the top quark decay. For illustration, we show 
in Table U the predictions in the presence oi a Wtb dipole 
interaction of the form 

- gRbLcr^’'tRdyW- , (18) 

whose strength is driven by the anomalous coupling gn. 
This amplitude is particularly interesting because it in¬ 
terferes with the standard model amplitude in the mass¬ 
less b limit. We select two values for this parameter, 
real and purely imaginary, which are well within the Icr 
limits given by current measurements (2^, [ 2 ^ . As it is 
apparent, the measurement of these additional polariza¬ 
tion observables will also improve the sensitivity to new 
physics in the top quark decay. 

Application to diboson resonances. For a heavy reso¬ 
nance R decaying into WZ 01 WW the W boson spin 
observables can also be measured using the method out¬ 
lined here. We consider as example a heavy scalar or 
vector resonance produced in quark-antiquark annihila¬ 
tion, focusing for definiteness on W'^ boson decays. We 
set our z axis in the direction of the W'^ momentum in 
the centre of mass frame. The remaining two axes are 
fixed by taking the x axis in the production plane, and 
requiring that the momentum of the diboson pair in the 
laboratory frame (used as a proxy for the direction of the 
initial quark momentum) has a positive x component. 

We give in Table |ll] the spin observables for a reso¬ 
nance with a mass of 2 TeV, as the diboson resonance 
peak found by the ATLAS Collaboration Q , and a mod¬ 
erate width of 50 GeV. The numerical calculations are 
performed with the Monte Carlo generator Protos [^ . 
It is known that for both heavy scalar or vector diboson 
resonances the resulting WjZ bosons have a big longi- 
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tudinal alignment of almost 100%, with {S 3 ) = 0 and 
(To) ~ —2fy/Q ~ —0.816, i.e. they are mostly produced 
with zero helicity. We observe in Table [TTl that the dif¬ 
ferences in (Tq) between scalar and vector resonances are 
indeed tiny and most probably unmeasurable. However, 
(Si) vanishes for a scalar while is of the order of 0.1 for 
a vector. This difference could be exploited to determine 
the spin of the resonance, provided it is detected in the 
W ^ ii' decay channel, as it should be if the ATLAS 
excess indeed corresponds to a. WZ or WW resonance. 

To summarize, in this Letter we have set a framework 
for the study of all spin observables of the W boson from 
its decay. We have given two examples of the usefulness 
of these observables, for W bosons from the decay of po¬ 
larized top quarks and for diboson resonances. The po¬ 
tential applications are numerous, as the formalism pre¬ 


sented here is completely general. Because W bosons are 
frequent decay products of (new) heavy particles, this 
formalism is expected to be of great value in the search 
of new physics effects. 
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